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Abstract. Two-term Weyl-type asymptotic law for the eigenvalues of one-dimensional 
fractional Laplace operator (— A)"/^ (a e (0,2)) in the interval (—1,1) is given: the n-th 
eigenvalue is equal to (n7r/2 — (2 — Q!)7r/8)" + 0{l/n). Simplicity of eigenvalues is proved 
for a G [1,2). and properties of eigenfunctions are studied. We also give precise 
numerical bounds for the first few eigenvalues. 



1. Introduction and statement of the result 

Let D = (—1, 1) and a G (0, 2). Below we study the asymptotic behavior of the eigenvalues 
of the following spectral problem: 

■—] <^{x) = X<^{x), xeD, (1) 

where G L'^{D) is extended to R by (for details, see below). It is known that there 
exist an infinite sequence of eigenvalues A„, < Ai < A2 < A3 < and the corresponding 
eigenfunctions form a complete orthonormal set in L^{D). The following is the main 
result of this article. 



Theorem 1. We have 



A„=(!!^-(^V + o(i), (2) 



2 8 y 

More precisely, there are absolute constants C, C such that 

C{2-a) I 



nix (2 — a)7r 



< 



a n 



forn> (C7«)3/(2°). 

The scaling property of the fractional Laplace operator {—cP/dx"^)^ implies that Xn{kD) = 
k''"'Xn{D). Hence, one easily finds the asymptotic formula for any interval. 

By following carefully the proof, one can take e.g. C = 30 000 and C" = 4 000 above. 
Note that the constant in the error term 0{l/n) tends to zero as a approaches 2, and in the 
limiting case a = 2 (not considered below), we have A„ = (n7r/2)2 without an error term. A 
stronger version of Theorem [1] for a = 1 was proved in |13| . 

The proof of Theorem [1] is modelled after [13]. In Section |2l an estimate for the fractional 
Laplace operator is given. The formula for the eigenfunctions on the half-line from [H] is 
recalled and studied in Section [21 An approximation to eigenfunctions is given in Section HJ 
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\2 


A3 
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u.yyi 


1.009 


1.009 


1.014 


1.014 


U.i 


u.yoi 


u.y/o 


1 riQi 
i.uyi 




1 1/17 


1 1 /IS 




Q71 


u.iyo 1 


1.195 


1.197 


1.319 


1.320 


0.5 


0.991 


0.970 


1.598 


1.601 


2.029 


2.031 


1 


1.178 


1.158 


2.749 


2.754 


4.316 


4.320 


1.5 


1.611 


1.597 


5.055 


5.059 


9.592 


9.597 


1.8 


2.056 


2.048 


7.500 


7.501 


15.795 


15.801 


1.9 


2.248 


2.243 


8.594 


8.593 


18.710 


18.718 


1.99 


2.444 


2.442 


9.733 


9.729 


21.820 


21.829 



Table 1. Comparison of the ap- 
proximation ~Xn ^ ^2;^)" 
(roman font), and numerical ap- 
proximations to Xn obtained using 
the method of [T5] with 5000 x 5000 
matrices (slanted font). 



Theorem [T] is proved in Section [5l and three further properties of eigenfunctions and eigenval- 
ues are studied in Section [6l Sections |3H6] correspond to Sections 8-10 in |13] . Proposition [3] 
gives the simphcity of the eigenvalues when a G (1,2). The result follows relatively eas- 
ily from the result for a = 1 in [13j. In Propositions [1] and [21 L'^{D) and L°°{D) bounds 
for eigenfunctions are given. Finally, in Section [TJ numerical estimates of A„ in terms of 
eigenvalues of large dense matrices are obtained. 

First-term Weyl-type asymptotic for A„ was proved by Blumenthal and Getoor in 1959 [3] . 
The best known general estimate for A„ is |(^)" < A„ < (^)" due to DeBlassie [9] and 
Chen and Song [7]. The important case of a = 1 was studied in detail by several authors, 
see [T| IT3] and the references therein. It is known that (A^i)^'^" is continuous and increasing 
in a G (0,2], see [21 [HI O [10]. For a discussion of related results and historical remarks, 
see e.g. [H [T3]. Theorem [T] is of interest in physics, the asymptotic formula (without 
the information about the order of the error term) was supported by numerical experiments 
in [TB], and there is a considerable amount of related (mostly numerical) research in physics 
literature. 

Noteworthy, although the values of C and C given above are rather large, numerical 
evidence suggests that the error term in formula ([2]) is rather small also for small n in the 
full range of a G (0, 2), see Table [1] and the estimates in the last section of this article. It is an 
interesting open problem to prove Theorem [1] with C and C non-exploding as a approaches 
0. This is related to simplicity of eigenvalues A„, conjectured to hold for all a G (0, 2), proved 
for a = 1 in [13j, and extended to a G [1,2) in Proposition [3] in Section [HI 

Motivated by the results of [TS] and [2], as well as by Theorem [T] above, one can conjec- 
ture asymptotic law similar to ([2]) for eigenvalues on an interval for more general operators 
A = ip{—d'^/dx'^), studied in |J^. While such a result for each individual ip should present 
no difficulty (under some reasonable assumptions on the growth of ip at infinity), it is an 
interesting (and much more dificult) problem to obtain estimates uniform also in ip, for a 
given class of ip. One important example here is the family of Klein-Gordon square-root 
operators A = — cP/dx'^ — m, with mass m ranging from to oo. This operator is 

close to \/ —d^ jdx^ for small m, but when m is large, it more similar to —(^jdx^. 
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To give a formal statement of the spectral problem ([T]), we recall the definition of the 
one-dimensional fractional Laplace operator A = {—d?/dx'^)"/'^. It is defined pointwise by 
the principal value integral, if convergent, 

Af{x) = c„ pv r ^^"^^" {j^^ dy, xeR, (3) 

^— oo p y\ 

where 

Af{x) is convergent if, for example, / is smooth in a neighborhood of x and bounded on R. 
Note that 

|a(2 -a) < Ca < |a(2 -a). (4) 

For / G C^(R), the Fourier transform of Af is equal to ^ extends to an 

unbounded self-adjoint operator on L^(R). We write Ad for the operator A on D with zero 
exterior condition on R \ More precisely, for / G C^{D), AdI is defined to be the 
restriction of Af to D. Again, Ad extends to an unbounded self-adjoint operator on L'^{D). 

The operator —A (on an appropriate domain) is the generator of the one-dimensional 
symmetric a-stable process Xt, and —Ad is the generator of Xt killed upon leaving the 
interval D. This probabilistic interpretation is a primary source of our motivation, but will 
not be exploited in the sequel. 

Notation. Throughout this article, C denotes an absolute constant (independent of a). We 
will track the dependence of other constants employed below on a to catch their asymptotic 
behavior as a \ and a 2. For brevity, we denote P = 2 — a. 



2. Auxiliary estimates 
Define, as in [13j, Appendix C, an auxiliary function: 



q{x) 



lix + l)^ 
-t(x 



1)^ 



for X G (— oo, 

for X G (-|,0), 

for X G (0, |), 

for X G (|, oo). 



(5) 



Note that q is piecewise C^, and q{x) + q{—x) = 1. Fix a a piecewise function / on R, 
and let g{x) = q{x)f{x). Further, we assume that the support of g is compact. Below we 
estimate Ag on (—1,0) in a very similar way as in [13]. 

Choose M to be the supremum of max(|/(x)|, |/'(x)|, |/"(x)|) over x G (— |, |). Let 
/ = \f{x)\dx. Then 



\g"{x)\ < \f{x)q"{x)\ + 2 |/'(x)g'(x)| + \f"ix)q{x)\ < CM. 
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Suppose first that x G (— 1, — |). Since g vanishes in (— 1, — i), c^^|^5f(x)| is bounded above 

by 



\9{y)\ 



1 \x — y 



l + Q 



dy <M 



q{y) 



)l+a roo 



< 



2i-"3"Mo 3 



2 - a 



+ 



dy 

l+aj 



21+c 



< 



CM 



\f{y)\dy 



/3 



+ CI. 



In the second inequahty we used the estimate — ^l^"*"" < + For x G (— 1,0) 

the principal value integral in the definition of A can be estimated by splitting it into two 
parts. By Taylor's expansion of g, we have 



pv 



gjx) - g{y) 



dy 



< sup (I \g"{y)\ : y G {x-\,x + \)] 
^ ^nv{\g"{y)\ : |/G (-|,|)} ^ CM 



'3 {x-yf 



\x — y\^^°' 



dy 



32-"(2-a) - 13 

Here for the second inequality note that g"{y) = for y < — |. Furthermore, 

9{y) 



We conclude that 




dy + 3 



l+a 



\f{y)\dy<^ + CI. 



a 



c-'\Ag{x)\< 



X G (-1,0). 



(6) 



3. Estimates for half-line 

The main result of [14j is the formula for generalized eigenfunctions for a class of operators 
on (0, oo). The case of fractional Laplace operator is studied in [14J, Example 1. In particular, 
the eigenfunction Fx of ^(o.oo) corresponding to the eigenvalue A" (A > 0) is shown to be 
Fx{x) = F{Xx) = sin(Ax + ^) — G{Xx) (recall that /3 = 2 — a), where G is a completely 
monotone function. More precisely, G is the Laplace transform of 



7(s) 



2a sin(^f ) 



27r 



l + s2°-2s"cos(^f) 



exp 



1 + 



log- 



1^ 



dr 



Furthermore, by |14| . Lemma 13, we have 

G{s) < sin(^) < C/3 

and 

/■oo 



poo 

/ G{s)ds = cos( 
Jo 



< C(5. 



Note that the exponent in ([7]) is negative. Furthermore, for a G (0, 1] we have 



(7) 

(8) 
(9) 



1 + s"" - 2s" cos(^f ) > (sin(f^))' > a% 
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while for a G (1, 2), the left hand side is not less than one. Hence, for all a G (0, 2], 

1 + s^" - 2s" cos(^) > min(a^ 1) > ^ • 
Finally, sin(^) <a(2 — a) = a/3. Therefore, 

2y2a/3 „ , , 

l{s)<— -s". 10 

avr 

By direct integration of the Laplace transform, we obtain that 



an \ a 



In a similar manner, (ITU]) gives 



< G"{s) < ^ (12) 

\/a Ja 



4. Approximation to eigenfunctions 

Let n be a fixed positive integer and /i„ = ^ — Our goal is to show that /^^ is close 
to A„. Note that /i„ > f and ^ < /i„ < ^. 

We construct approximations to eigenfunctions by combining shifted eigenfunctions 
for half-line, F^^(l + x) and -F)j„(l — x), and using the auxiliary function q given above in ([5]) 
to join them in a sufficiently smooth way. We let 

^„(a;) = g(-x)F^„(l + x) + (-l)"g(x)F^„(l - x). (13) 

Lemma 1. We have 

cp> 1 

||^D<^n-/in<^n||2 < (14) 

A/a n 

Proof. Note that we have 

^„(x) - F^„(l + x) = -(1 - g(-a;))F^„(l + x) - (-l)"g(x)F^„(l - x) 
= -g(x)(F^„(l + x) + (-1)"F^„(1 - x)) 

= g(x)(G^„(l + x) + (-1)"G^„(1 - x)) - sin(/x„(l + x) + f )l[i,oo)(x). 



X 



Denote ;i(x) = sin(/i„(l + x) + ^)l[i,oo)(x) and /(x) = G^„(l + x) + (-1)"G^„(1 
(7(x) = g(x)/(x). It follows that (^n(a;) = F^„(l + x) + g{x) + h{x). For x G (—1,0), we have 
AFfj_^{x) - /i"F^„(x) = and h{x) = 0. Hence, 

\A^n{x) - fi':^n{x)\ < \Ag{x)\ + \Ah{x)\ + IfiMl X e (-1,0). (15) 

We will now estimate each of the summands on the right hand side. 
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Using convexity of G, — G" and G", and estimates ([9]), (fTTj) and ( IT2|) . we obtain that 
sup |/(a;)| < + < ^/x;^"", 



sup < -^^G'd/i.) -/inG"(|/i„) < ^/i^^-", 

sup |r(x)| < < ^/i;;^-" 



oo />oo 



\f{x)\dx< I Gi_,^{l + x)dx + J Gi^^{l-x)dx 
^ rG{y)dy<^. 

f^n Jo f^n 



By dn]) and 



Furthermore, \g{x)\ = for x E (—1, —\), and 



^^(x) I < ^ + X e (-1, 0). (16) 



K9{x)\ < yUix) < X e (-i 0). (17) 



'a 

Finally, for x < we have the following estimate for the oscillatory integral 



\Ah{x) 



sin(/i„(l + |/)+ (g ^- 
\x - 



8 



dy 



- Ji |x-l|i+° ^ - (1 - x)i+"/i„ - "^^^ 

Estimates ([I6D-([I8]) applied to yield that 

CB 

\Aipn{z) - fil^n{z) I < 2; e (-1, 0). (19) 



By symmetry, ( fT9l) also holds for 2; G (0,1). Formula dH]), with ^D<^n understood in the 
pointwise sense, follows. It remains to prove that (pn is in the domain of Ad- To this end, 
we will use the notion of the Green operator Gd = A]j^. The reader is referred e.g. to [6] 
for formal definition and properties of G/j. 

Since Aipn is bounded on D, the function <^„ — GDA(fn is a bounded, continuous in D, 
weakly a-harmonic function in D = (—1, 1) with zero exterior condition. Such a function is 
necessarily zero (see [U [H]). It follows that = GoAifn, and hence (fn is in the L°°{D) 
domain of Ad- Since convergence in L°°{D) is stronger than the one in L'^{D), the proof is 
complete. □ 



Lemma 2. We have 



CB „ „ Cf3 

— < <^n 2<1 + — . 20 

n n 



In particular, there is an absolute constant K such that \\'^n\\2 > \ for n> K . 



~ /in 



EIGENVALUES OF THE FRACTIONAL LAPLACE OPERATOR IN THE INTERVAL 
Proof. First, note that by direct integration, 

(sin(/i„(x + + - dx 

Using fll3p and we obtain the lower bound, 

> j {sm{fin{x + 1) + f)Y dx 

+ l)sin(/i„(x + 1) + f)| c/x > 1 . 

-1 * ' /in 



I ~ Il2 



In a similar manner, 



I ~ Il2 



< 1 + — + 4 [ {Gifinix + < 1 + ^ 

/in J-1 /in 



and the lemma is proved. □ 

5. Proof of Theorem [T] 

Since ipn € L'^{D), we have ipn = Ylj^^jVj ^'^^ some aj. Moreover, ||<^n||2 = '^j ^'^'^ 
^D'Pn = X]j ^j^j^j- Let Afc(„) be the eigenvalue nearest to /i". Then 



oo 

a\2\\,~ ||2 
n||2- 



i=l J=l 

By (I14p and Lemma HI it follows that for n > K, 

\>^Hn)-l^:\<^-. (21) 



This will enable us to derive a two-term asymptotic formula for \j. 
Denote £ = ^ We have 

|(/i„ ± - /i^l > aemin((/i„ - + e)"-^) > Caen"-\ (22) 

Thus if |Afc(„) - /i^l < then An G ((//„ - e)", (/Xn + e)")- By dH]), this holds true if 

n > K and 

Therefore, Lq, = |"(Ca~^/^)^/"] (the constant here is chosen so that also La > K) is such 
that for n > L^, each interval ((/i„ — e)", (/!„ + £:)") contains an eigenvalue Xk{n)- In particular 
Xk{n) are distinct for n> L^. We claim that there are less than eigenvalues not included 
in the above class. As in [T3], the key step will be the trace estimate. 

Let J be the set of those j > for which j ^ k{n) for all n > La- Denote by pt{x — y) 
and p^{x,y) the heat kernels for A and Ad respectively; we have pt{C,) = exp(— 1|^|"). For 
t > 0, we have (see e.g. P 

Ve-^^*= / Ve-^^*((^j(x))2c/x = / p^{x,x)dx< / pt(0)dx = - / e~*^°ds. 
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In the last step, Fourier inversion formula was used. Hence, 

jeJ j=l n=La ^ n=La J 

The latter series is bounded below by the integral of e~*'^° over {fiia + ^, oo). Hence, 



2 rf^La+s 2 



Taking the limit as t \ 0, we obtain that 

„ . 2 , , ^ (3 2e 

#J< - (/XL. +£) = ^a-7 + — • 

TT 4 77 

Since e < the right hand side is less than L^, and the claim is proved. 

By [HI |9], we have A„ < (n7r/2)". It follows that for all n < La, we have A„ < (/i^^ — e)", 
and so J = {1,2, ...,Lq, — 1}. We conclude that k{n) = n for all n > La- Theorem [1] now 
follows from ( 12T]) . 

6. Further properties of eigenvalues and eigenfunctions 

Int this section three additional properties of (fn and A„ are studied. This part is modelled 
after pj3j. Section 10. A number of open problems is suggested at the end of the section. 

Proposition 1 (cf. Lemma 3 and Corollary 4 in [13j). There is a constant C such that, 

C{2-a) 

\Wn-'^nh< whena>l, 

n 

„ . „ ^ C{2-a) 

\Wn-Vn\\2< — 0/9 a whena<l. 

In particular, if ip*^{x) = ±cos(/i„x) for odd n an fni^) — ±sin(/x„x) for even n, then 

> * u C{2-a) , 1 

Wn-^nh< 1= — whena>-, 

In I 



,. C(2-a) ^ 1 

- ^nh < ^3/2^a ^^^^ " < 2 ' 

for some constant C . 

Proof. Fix n > La + 1 and ^ = \ and write, as in the previous section, (pn = Ylij ^j^j- 
Changing the sign of ipn if necessary, we may assume that a„ > 0. As in ( 122|) . for j ^ n we 
have \\j — > Carf~^ . Hence, 

oo 

j=l jjtn 

By (fT^ . we obtain that 

jj^n \ V / \ / \ / 
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Note that 

WVn - \\Vn\\2Vn\\2 < WVn " «nV^n|l2 + " ll'^nlhl • 

But |a„ — ||<^n||2p < ll'^nlli ~ '^n = ll*!^" ~ Qn'/'nlli- Heiice, using also (120|) . we obtain that 

2C/3 , 

WVn-WnhVnh ^ ^zj^ ' (24) 

Finally, by ([23]), 

,. ., , 2C(5 CP 

We have thus proved the first part of the proposition. The second statement is a simple 
consequence of the first one, the identity (pn{x) = v^* (x) — G{fin{l + x)) ± G'(/i„(l — x)), and 
the estimate \\G\\l < ||G||i||G||oo < □ 



Proposition 2 (cf. Corollary 5 in [13]). If a > \, then the eigenf unctions ^Pnix) are bounded 
uniformly in n > 1 and x E D. 

Proof. Let P/^ = exp^—tAo) {t > 0) be the heat semigroup for Ad (or transition semigroup 
of the symmetric a-stable process in D), and let p^{x,y) be the corresponding heat kernel 
(or transition density). It is well known that p^{x,y) < Pt{y — x), where pt{x) is the heat 
kernel for A and pt{C) = 6xp(— see e.g. [5]. 

By Cauchy-Schwarz inequality and Plancherel's theorem, we obtain 



h /"e"'*'''"^'^ll¥'--<^"ll2+ Halloo 

< 2v/r(l + l/a) (2t)-V(2-)||^„ _ + 2, 



Let t = 1/A„. Then e"^"* = 1/e and t-V(2") < ^1/(2^) ^ If n > L« and a > |, then 

also \\ipn ~ 'Pn\\2 ^ C'/'y/n, and finally |<y9„(x)| < C. Since each is in L°°{D), the proof is 
complete. □ 

Proposition 3 (cf. Theorem 6 in |13j). //«>!, then the eigenvalues A„ are simple. 

Proof. Let us write A„,q, for in this proof. Since (A„,q,)^/" is increasing in a, we have 

(A„,„)^/" < (A„,2)'/' 



T ■ 

By Theorem 6 in [13], for > 3 we have 

^ 8 ~ 10 - '^^"+i.aJ 

Therefore, A„_q, < A„+i.q,, except perhaps n = 1 or n = 2. But a similar argument works also 
for n = 1 and n = 2, since by [T] we have 

TT 

- < 2 < A2,i, and vr < 3.83 < As^i. 

2 



10 
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The proof is complete. □ 

Numerical experiments suggest that are uniformly bounded also for a < |. Further- 
more, it would be interesting to obtain an upper estimate of sup„ llv^nlloo, and in particular, 
to find its behavior when a approaches 0. Finally, as stated in the introduction, better 
bounds for A„ may yield simplicity of eigenvalues also when a < 1. 



7. Numerical bounds for eigenvalues 

No general efficient algorithm giving mathematically correct numerical bounds for A„ is 
known to the author. For a = 1, a satisfactory method (an application of Rayleigh-Ritz 
and Weinstein-Aronszajn methods) is described in [13j. For general a, even approximation 
of An is difficult: all known methods converge rather slowly, and thus the computation of 
eigenvalues of very large matrices is required. In this section a method for obtaining a 
lower bound for A„ is described. It shares the main drawbacks of many related algorithms: 
compared to the technique applied in [15j, it converges slowly, and it suffers large errors 
as a approaches 2. On the other hand, the method presented below gives mathematically 
correct lower bounds, and there is no error estimate for the numerical scheme of [15]. At the 
end of the section, a somewhat similar method for the upper bound for Ai is given. It gives 
satisfactory results for large a, but deteriorates as a gets close to 0. 

It should be pointed out that in many particular cases [a close to 2 or large), the bound 
|(^)" < A„ < (^)" of [71 [H] is sharper than the estimates obtained below, unless extremely 
large matrices are used. Also, good numerical estimates of A„ are available for a = 1 due 
to [13], By the monotonicity of (A„)^''" in a, this gives a lower bound for A„ when a G (1, 2) 
and an upper bound for a G (0, 1). Finally, a good estimate of Ai can be found in [1]. For a 
comparison of the above, see Table |2l 

Our method for the lower bound works for fractional Laplace operator in an arbitrary 
bounded open set D C R"^ (in fact, it can be easily extended to more general pseudo- 
differential operators, or Levy processes). Fix e > and let |Jfc : A; G Z"^} be the partition 

of R"^ into cubes h = 11^=1 i^j + 1)^], k G Z^. Let C Z"' be the set of those k e Z'^ 
for which intersects D, and let D,, be the interior of IJfce-K' ^'^^^ that D <Z D^. 

The definition of ^ = (— A)"/^ in higher dimension is similar to ([3]): for smooth bounded 
functions we have 

\x-y\'^+'^ 

where Cd^a = 2"r(((i + a)/2)/(7r'^/^|r(— 1)|). Fractional Laplace operator in D with zero 
exterior condition, denoted Ad, is defined as in dimension one. Below we denote by A„ the 
eigenvalues oi Ad- By domain monotonicity of A„, the eigenvalues for D are not less than 
than the eigenvalues of its superset D^. For notational convenience, we assume that D = D^. 
The Dirichlet form S{f, /) corresponding to Ad is given by 
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A, 




'^3 






'^6 




^8 


An 






u.yyoD 


i.UUoD 


i.Uio ^ 


1 ni Ti 
i.Ui ( i 


1 ni Afi 


1 A01 7 
i.Uzi / 


1 noQ /i 


1 no/1 Q 


1 no^i 1 
i.UzDi 


1 AO70 




u. yy4o 


U.OUO / 


U.OU/o 


u.ouyz 


A fr 7 A/f 2 
U.01U4 


A c: 7 7 Q-^ 
U.Ollo 


A f: 7 07 
U.OIZI 


A fr 7 ooS 
U.OlZo 


A fr 7 Q/f 2 
U.01o4 


u.oioy 


0.01 


(J.997D 


1.0086 


1.0138 


1.0172 


1.0198 


1.0218 


1.0235 


1.0250 


1.0263 


1.0274 




l.UUo/ 


7 ni Q74 
l.Ulo/ 


7 D 7 704 

l.ul iZ* 


7 D 7 7^ 

l.UW / 


7 noi q4 
l.UZlo 


l.UZoO 


l.UZOU 


l.UZOo 


7 /)07A4 
1 . UZ / 4 




1 Q KOI n 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 




0.9974^ 


1.0102'^ 


1.0148^ 


7 ni7Q^ 


7 090'^''^ 


i 0223^ 


7 n9'?q'^ 


1.0254'^ 


7 0966''^ 


7 0977^ 




U.y / z4 


i.uyiy 


11/1 ^^o 


1 1 SfiQ 
i.loDo 


1 01 p;o 
i.Zioy 


1 0/inr^ 


1 9(^1 1 
i.ZDii 


1 0701 

i.z / yi 


1 90p;n 

i.zyou 


1 Qno/1 




(7. yoio 


U.OOUO 




U.OUUo 


A f;-} AA^ 
U.0144 


U. OzD / 


U.O0O4 


U.044U 


A /?fr 7 


U.OOoO 


0.1 


r\ nonn 


i.uyio 


11/1 TT 

1.14^^ 


1 1 OCT 


1 01 i^T 

l.ZlD / 


1 /1 1 
l.z41Z 


1 oc;oA 
l.ZDZU 


1 OOAO 

l.zoUz 


1 0Ai?0 

i.zyoz 


1 01 AT 


U.hf / ZO 


i . UaZZ 


7 7 7'^4 

1.14:/ o 


7 7 S^s4 
l.lOOO 


1 .ZlOO 


7 041 '^4 
1.Z410 


1 .ZOZU 


1 .zouz 


7 00(^o4 
1 .Z\yOZ 


7 1/174 
l.olU / 




l.OOOi 


n/a 


n/a 


n/ a 


n/a 


n/a 


n/ a 


n/ a 


n/ a 


n/a 




0.9786'^ 


1.1067'^ 


1.1575^ 


i i 941 ^ 


1.2226^ 


1.2462^ 


1.2664^ 


1.2840^ 


1.2997^ 


1.3138^ 




U.yo / z 


i.iyou 


1 Q 1 Q9 
i.oioZ 


1 /inoQ 


1 /I QOI 
i.4oUi 


1 c:/ino 


1 p;oi p; 
i.oyio 


1 fiQ7Q 
i.DO / 


i.D / oU 


1 71 Ci/I 
1. ( i04 




u.yioi 


U.ozoO 


U.Ool / 


A TOO 7 S 


U. /OOU 


A TO Q 7 

U. /ool 


U.oU/0 


A QOOA^ 

u.ozy4 


A 0/7A02 

u.o4yz 


U.oO/o 


0.2 


n ATI o 

U.y iz 


1 1 n /I o 

i.iy4o 


1 Q 1 nn 

i.oiyy 


1 /II AO 

1.41Uz 


1 /I 1 A 

i.4oiy 


1 K /I OA 

1.04ZU 


1 K AQA 

I.oyoy 


1 ^;OAA 

i.Doyy 


1 ^;o 1 
l.DolZ 


1 T1 

1. 188 


u. yo/0 


7 7 Qft^4 


7 Q1 Ql4 

i.oiyi 


7 Ain^4 
1.41U0 


7 AQ1 74 
1.4ol / 


7 KA0'\4 
1.04Z1 




1.04UU 


7 7 7-4 

l.Ooll 


7 71 qq4 
1. / loo 




i.zo / D 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 






1.2247^ 


7 '?*?Q«'5 
± .OOc/O 


1.4258^ 


1.4947^ 


7 

J. .'JOOU 


J- .uuou 


7 fid^^i^ 

1 .U'rtOO 


7 fi^QO^ 

1 .\jOc/U 


7 79(^0^ 




U.yDyz 


i.oyyi 


o no /1 7 
z.Uz4 / 


QOAO 


Z.DODZ 


A^il Q 

z.yoio 


Q 01 1 Q 

o.ziio 


A A AQ 

o.444o 


fi^iAQ 

o.dOUo 


Qt; cr /I 
O.OD04 




U.OoOZ 


U.OoOZ 


7 rtQi^4^ 
1 . Uo04 


l.ZOoo 


7 ACilO^ 
1.4U1Z 


i.0o4y 


1.00/ y 


7 779//^ 

i . / / z 4 


1.0 /yy 


7 OQIfi^ 

i.yoio 


0.5 


u.yyuo 


i.oy / / 


Z.UoUd 


ooco 

Z.OODZ 


^JACC /I 

z.Dyo4 


ATO K 

z.y /^Zo 


OOKA 

o.zzoy 


/I c; AO 


f^OAO 


0000 
0.0000 


U.y/Ui 


l.OUlo 


o noQQ4 
z. UZoo 


Z.oo/1 


/^A//7^ 

z.oy4 / 


z.y/zo 


o.zzOO 


/1 1^ 7 a4 


o.ooUO 


ooco-^ 
0.0000 




l.UUUz 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 




0.9863^ 


1.6598^ 


2.0777^ 


2.4274^ 


2.7314^ 


0. uuoo 


3.2562^ 


3.4892^ 


3.7074^ 


3.9136^ 




l.lOiD 


Z. / o4o 


A '~>7^fi 
4.Z / OD 


O.oZoO 


7 qp;q/i 
/ .OOo4 


Q Q AI A 


1 n /1 1 ^ifi 


1 1 OQQO 


■\ A C^OQ 
iO.40Zo 


1 /I OfiQfi 
i4.yDOD 




1.10/ / 


9 7^4 7^ 
Z. / 04 / 


4.0100 


cr 0097 5 
O.ObfZl 


7 Aftm 

/ .40Ui 




7 A fifiOO^ 
1 U. OUZZ 


7 9 1 741 3 
IZ.l / 41 


70 7441 3 
lo. / 441 


7 fr 7 f; c: 5 




1.1781 


o ^ /I on 

Z.7489 


4.3197 


5.8905 


7.4613 


9.0321 


1 A i^AOA 

1U.dUz9 


101 TO T 

lz.1737 


13.7445 


15.3153 


1 


7 7 K77-^ 
1.10/ / 


9 7^4 ^4 
Z. / O'lO 


4.0 104 


O.ObflO 


7 4^044 
/ .40bf4 


A n 9 7 


7 n ftn 7 9^ 

lU.OUlZ 


797 790-^ 

IZ.l / zy 


7 74074 
lo. / 4Z / 


7 fr 07 JA^ 
10.O14U 




l.lOUo 


n/ a 


n/a 


n/ a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 




7 7 cjrfi-^ 


9 7^dR^ 




O.Oc/^^ 


7 4609^ 


Q 0'^9Q^ 


J- L/. UU^O 


12.1742^ 


13. 7442^ 


7 Q7 c;^,'^ 




i.oioy 


4. / oD / 


Q QQ1 7- 


iO. / DOO 


1 ot^no 
ly.zouz 


OK 0(^1 Q 


01 7QQ/1 

ol. / oo4 


00 fJOf^Q 


4o.oyyD 


fro c:o(^^^ 
Oo.OZdD 




1 . OZiJo 


A 79 7 '■^ 
4.0/ zl 




14.0UZ4 


OA Q7fiO°^ 

ZU. / oz 


97 7 A70^ 

z/.i4/y 


Q/i r.ooo^ 
O4.0ZZZ 


/fo A 770°^ 
4Z.4/ / Z 


fTA QCiOfi^ 

OU. yooo 


fTA AQTC;5 

oy.yo/o 


1.5 


1 C 1 1 /I 

1.0114 


o.Uo4o 


A C ATA 

y.oy /U 


1 K A1 TI 
lO.Ul i 1 


01 1 AAC 

zi.iyuo 


00 AO /I /I 

zo.Uo44 


OK /I OOC 


A c: ac;t 


CO AK 1 A 
OZ.UOI4 


CI AAOO 

Di.uyzz 


l.OiJ / 1 


O.UOoO 


y.OiJzi 


1 p: m ka4 
10.U104 


97 1QA/^4 
Zl. Io40 


00 noQo4 
Zo.UZoiJ 


o0.4oUU 


A0704 
4o.4y / Z 


fro nQQ04 

OZ. uoyz 


^?7 r\7Qfi4 

Ol.U/oO 




i.oyoy 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 




1.6224^ 


'J.tJUO~± 


7 9907^ 


7 c 7dQ7^ 
10.1 / 


99 mn^^ 

iLi^. Ij-L iJO 


9>? O'^'^q^ 


*?/^ 4600^ 


izl ^id67^ 


^'^ 7 crcrn^ 

J- !j<Jl7 


/^9 9K^fC^ 




1.44oo 


1 1 /IQ 

0. Il4y 


in AAA7 
1U.444 ( 


1 T 00Q1 


OK 0Qn7 


OA Pi/1 /IS 
O4.044o 


A A SQf^Q 

44.oyDy 


K« 9S1 Q 
OD.ZolO 


CO fiOQc; 
Oo.DooO 


01 Qoin 
oi.yziu 




1.0 / 00 




lo.yuoo 


OA QAOft^ 
Z4.o4bfO 


07 0QA7^ 
o/.Zo4/ 


cro croAQ^ 


TA 7 AAo5 


OA AAr7^' 

ou.yuo / 


7 7 7 SAQO'^ 

111. o4oz 


7 Qfr Cir\fir\3 

loO.yuou 


1.8 


2.0555 


7.5003 


15.8014 


0(? TO 00 

ZD. 7233 


4U.1148 


55.8658 


TO OAAC; 

73.8905 


A /I 1100 

94.1188 


11c -1 AO 

110.4923 


1 /I A AC A f; 

140.9605 


Z.Wlol 


/ .ouu/ 


7 ^ 7QAS4 

10. / y4o 


ZO. / 100 


A 1 m o4 

4U. lUlZ 


crc: Q4Qi4 
00.O4O1 


70 Q(ifil4 
/ o.oOOl 


A/7 nQQA4 
iJ4.U004 


77c 4 1^41 4 
110.4041 


7 /7n 01 4^4 

14U.y 140 




Z.UOUi 


n/ a 


n/a 


n/ a 


n/a 


n/a 


n/ a 


n/ a 


n/ a 


n/a 




9 0777^ 

jC.U / / / 


7 


±u.^ouo 


97 '^'^^'^^ 


z/n >?J79^ 
■J I/. 0^ / 


W 77 "^^^ 

OU. / ±00 


74 R^iOl^ 


Qc: 1R71^ 
i/0. 1 i ± 


7 7 7 fifi/^4^ 


142.2381^ 




l.Uooo 


6. i / U4 


7 Q7Q/1 


101 OQO 


1 (^Q70 

ly.Do / y 


07 1 1 KG 

z / .lloy 


oo.ooyi 


A A 0/1 CI 
44.y4oi 


CC OHQO 

OO.ZUoZ 


cc q 1 OT 
DO.OlZ / 




7 C97Q^ 
i.OZ/O 


« Q^7Q'^ 
O.oO (o 


7 /^OOQ-^ 

i 0. uyuo 


OA n7Kn^ 
Zif.U/OU 


/7 ^001 °^ 
40. OZZl 


^c; A7Q7^ 
00.4/0/ 


00 7(iQfi3 

00. /OoO 


7 7 r AQQQ'^ 
110.4OOO 


7 fr Q c; 7 5 
140.O0Z1 


7 7C KAfiQ^^ 
1 /o.040o 


1.9 


2.2477 


8.5942 


10 Tl TT 

18.7177 


32.4615 


AC\ TO A /I 

4y.7zU4 


TA /1 1 T 

7U.4157 


A /I A A 

94.4848 


101 TC /I 

121.8754 


1 cr fr /I 

152.5433 


10c /If; AA 

18o.45U(J 


o 0/1 Qo4 


o.ObfZO 


70 71 m 4 
lo. / lUl 


00 A^nQ4 
0Z.4OU0 


AO 7noi4 
4y. /uzi 


TA QQne:4 


OA A'^nQ4 
iJ4.40Uo 


797 QQ1 q4 
IZl.oolO 


7 fC9 AQ7Q4 
10Z.4o/o 


7 QS004 
loO.ooZZ 




O 0/1 c: 
Z.z400 


n/ a 


n/a 


n/ a 


n/a 


n/a 


n/ a 


n/ a 


n/ a 


n/a 




2.2748^ 


o.ou^ ± 


7 q n 7 

J. i/. ui 1 


Q9 s^n^^ 


^rt 1Qfi9^ 

OU. 1 c/Uj^ 


7n Q766^ 


190'^^ 


799 fin9d^ 


7 c;q Qt:7 7^ 


7C7 '?'?«Q^ 
J-O / .000^7 




0.1474 


0.5494 


1.1671 


1.9816 


2.9788 


4.1482 


5.4811 


6.9705 


8.6101 


10.3944 




1.9816^ 


7.5 121-^ 


18.3642^ 


34.1070'^ 


54.5469'^ 


79.8163^ 


109.7856^ 


144.5508'^ 


184.0144'^ 


228.2517'^ 


1.99 


2.4441 


9.7330 


21.8288 


38.7113 


60.3666 


86.7839 


117.9546 


153.8713 


194.5275 


239.9178 


2.4427^ 


9.7293^ 


21.8200^ 


38.6960^ 


60.3426^ 


86.7495^ 


117.9077^ 


153.8100^ 


194.4500* 


239.8220-^ 




2.4452 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 




2.4563^ 


9.7573^ 


21.8651^ 


38.7595^ 


60.4267^ 




i 18. 0385^ 


153.9670^ 


194.6351^ 


240.0373^ 


^ See 


[IJ. ^See[7]. ^ 


Combination of 


"3] with monotonicity in a. 


^See [15]. 







Table 2. Comparison of bounds and approximations to A„. Each cell contains six 
numbers: lower bound A„^e with e = the best lower bound known before, ap- 



proximation (-^ — g"'*'^ )", numerical approximation of |15| . upper bound A* the 
best upper bound known before. The better estimates are printed in color. 
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As usual, / e L'^{D) is extended to R*^ so that f{x) = for a; e R*^ \ £). For k e Z'^, denote 

When X e Ik, y ^ lu k, I e Z*^, we have \x — y\ < £\\k — l\\. We define 
and 

Clearly, £e{f-if) < £{f,f)- By Rayleigh-Ritz variational principle, the eigenvalues A„ are 
bounded below by the sequence Xn,e oi eigenvalues of the operator corresponding to the 
Dirichlet form S^. More precisely, A„ ^ are defined in the usual way, 

Xn,e = inf {sup f):feL, Wf^ = 1} : L < L^D), dimL = n} . 

Here 'L < L'^{D)^ means that L is a linear subspace of L'^{D). 

For / G L'^{D) and k G Z*^, let = s^'^ Jj^ f{x)dx, and define /* to be equal to on I^. 
Hence /* G L'^{D) is the orthogonal projection of / onto the space of functions constant on 
each 7fc, and Jj^f*{x)dx = Jj^f{x)dx. In particular, ||/||2 = ||/*||2 + ll/~/*ll2- Furthermore, 

f) = E '"^-^ I I ((/(^))' - 2/W/(y) + {f{y)f) dxdy 

Comparing this with a similar formula for /*, we obtain that 

S,{fJ)=S,{f\n + c,,^e-^iy\\f-r\\l 

This shows that the two orthogonal subsj:)accs, {/ G L^{D) : /* = 0} and {/ G L^{D) : 
f* — /}) invariant under the action of the operator corresponding to S,.. The former 
subspace is in fact its eigenspace, corresponding to the eigenvalue Cd,a^~°'. The latter one is 
finite-dimensional, and the action of in the basis of normalized indicators oi 1^, k E K,,, 
is given by the following matrix V: ii k he a bijection between {1, 2, \Ks\} and K^, then 

We conclude that the sequence A„,£ starts with those eigenvalues of the matrix V which 
are less than Cd,a^^°'i^, which are followed by the constant c^^q^ This gives the lower 
bound for the eigenvalues A„ for an arbitrary open bounded set D. Note that replacing z/ 
be a smaller number gives smaller lower bounds A„_£, hence the series defining i/ should be 
approximated from below. 
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When D = (—1, 1) C R and e = j^, then u = 2({1 + a) — 1, where ( is the Riemann zeta 
function. Furthermore, in this case V is a Toephtz matrix with the symbol 

2crf,„ ( ^ ^ cos{kx) \ _2cd,a f^r, , „^ [ l^ii+aie'"') 



2Cda ( ^ 1 r t"(e*-COSx) , 



1 + « 7o e^* — 2e* cos x + 1 

The right hand side is easily checked to be increasing in x G [0,7r], and so it attains its 
maximum for x = n. The symbol of V, and hence the eigenvalues of V, are therefore 
bounded above by 2cd,a£~"'iCi^ + a) — Lii+a(— 1)) = 2^~"cd,a£~"Ci^ + a) < Cd,aE~°'i^- It 
follows that all N eigenvalues of V are included in the sequence Xn,e- 

In higher dimensions, u can only be computed by approximating numerically a d- 
dimensional infinite series. We summarize the results of this section in the following two 
results. 

Proposition 4. Let D = (—1, 1), > and e = 2/N. Let V be a N x N Toeplitz matrix 
with entries 

c 1 

^P^^ = + ' p, g = 1, 2, AT, p ^ g; 

^ 2c.(C(l + «)-!) 

l/pp , jj ...,iV. 

Define A„^e to he the n-th smallest eigenvalue of V when n < N, and Xn,e = CaB~°'{2({l + 
a) — 1) otherwise. Then the eigenvalues A„ of Ad satisfy A„ > Xn,e- 

Proposition 5. Let D C R*^ be an open set in H'^, and let e > 0. Let Kf. he the set of those 
k eZ'^ for which D n ]\%i[kje, {k^ + 1 )e] is nonempty, and let n : {1,2, \K^\} — )■ K,, be 
the enumeration of elements of . Finally, let 



^ \\k\\-'^-'^, where 
fcezd 



d 



\ i=l 



Define a \Kir\ x li^^l matrix V with entries 

^v,, = Mp) - <q)\\-''-'' , P, g = 1, 2, \K,\, p^q- 

= ^ - , p=l,2,...,iV. 

Let Xn^e be the n-th smallest eigenvalue of V if n < N and this eigenvalue does not exceed 
Cd,a£~°'^, OLnd Xn,e = Cd,a£~°'i' othcrwisc. Then the eigenvalues A„ of Ad satisfy Xn > Xn^e- 

The lower bounds A^^^ for the interval D = (—1,1) are presented in Table [2] above. In 
higher dimensions, the complexity of computations increases dramatically. For example, a 
unit disk -8(0, 1) or a square [— 1, 1]2 with e = ^ require handling matrices larger than 
2000 X 2000. Some results for these two cases are given in Tables |3] and |H 

In principle, the upper bound is much more difficult. The above approach can be modified 
to give an upper bound for Ai whenever the Green function for D can be computed. For 
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Oi 


M 




\ /'TTR^ 


A2 


(Lti) 


A2 {^ti} 


1 

U. J- 


1 O'^ns 

-L .UOUO 




1.0462^ 


1 ORRO 
± .uoou 




1 .wool- 


0.2 


1.0506 


0.5472^ 


1.0946^ 


1.1691 


0.5865^ 


1.1731' 


0.5 


1.1587 


0.6266^ 


1.2534^ 


1.4908 


0.7452^ 


1.4905' 


1 


1.3844 


0.7853^ 


1.5708^ 


2.1807 


1.1107^ 


2.2215' 


1.5 


1.4135 


0.9843^ 


1.9688^ 


2.6029 


1.6554^ 


3.3110' 


1.8 


0.9167 


1.1271^ 


2.2544^ 


1.8164 


2.1033^ 


4.2068' 


1.9 


0.5427 


1.1792^ 


2.3585^ 


1.0984 


2.2781^ 


4.5563' 



^ See [7]. 



Table 3. Comparison of estimates of A„ for 
a square [—1,1]^. LB and UB mean lower 
bounds and upper bounds respectively. Es- 
timates of this section are given in roman 
font, best numerical estimates known before 
are typeset in slanted font. Better estimates 
are printed in color. 



a 


Ai 


(LB) 


Ai 


(UB) 


Aa 


(LB) 


\2 (UB) 


0.1 


1.0381 


1.0157' 


6.6198 


1.064V 


1.0953 




1.1609'-' 


0.2 


1.0655 


1.0396' 


3.8878 


1.1342' 


1.1849 


0.6541'^ 


1.3476'^ 


0.5 


1.1986 


1.1618' 


2.5081 


1.3943' 


1.5404 


0.9787^ 


2.1079^ 


1 


1.4734 


1.5707' 


2.7588 


2.0944' 


2.3201 


1.9158^ 


4.4429^ 


1.5 


1.5387 


2.3891' 


4.0668 


3.4131' 


2.8379 


3.7502^ 


9.3648^ 


1.8 


1.0087 


3.2210' 


5.5014 


4.7468' 


2.0045 


5.6114^ 


14.6487^ 


1.9 


0.5990 


3.5834' 


6.1369 


5.2974^ 


1.2165 


6.4182^ 


17.0045^ 



^Seep. ^See[7]. 



Table 4. Comparison of estimates of A„ for a unit 
disk. LB and UB mean lower bounds and upper 
bounds respectively. Estimates of this section are 
given in roman font, best numerical estimates known 
before are typeset in slanted font. Better estimates 
are printed in color. 



the fractional Laplace operator, this is the case when D is a ball. By a scaling property, it 
is enough to consider D = B{0, 1). 

Let Gd be the Green function of D, Gd{x, y) = p^{x, y)dt, where is the heat kernel 
for Ad (see the proof of Proposition [2]). The Green function is the kernel of A]j^. M. Riesz 
proved that 



r ( \ - ^^^^ / rl 



2"7r'^/2r(f)r(l + f)|a;-|/|'^' ' V2' 2' 2' \x - y 

Since the eigenvalues of A]^^ are A~^, we have 

^ = sup 1^ ^ Gd{x, y)f{x)f{y)dxdy : / G L\D), 
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Since Gj:){x,y) is nonnegative, we may restrict the supremum to nonnegative functions only. 
Hence, whenever g{x,y) < GD{x,y), we have 

Ai < (^snp J jix,y)fix)f{y)dxdy:f E L\D), \\fh = 1 

For k,l G Z'^, let g^^i be the infimum of Gd{u,v) over u E Ik and v E Ii. When x G /fc, 
y G we choose g{x,y) = gk^i- Hence, Ai is bounded above by AJ^, the reciprocal of the 
largest eigenvalue of the matrix U with entries Uij = £'^5'K(i),K(j)- 

The results for D = (—1, 1) C R and some values of a are given in Table |2l Estimates 
for the unit disk and the square [—1,1]^ are given in Tables [3] and HI Noteworthy, for the 
unit disk and £ = the estimate A^ ^ is worse than the one obtained in [1] using analytical 
methods. 

Acknowledgments. I would like to thank Krzysztof Bogdan and Tadeusz Kulczycki for 
helpful discussion and valuable suggestions. 
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